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Abstract 

We perform the bifurcation analysis of the Kepler problem on and 
L^. An analogue of the Delaunay variables is introduced. We investigate the 
motion of a point mass in the field of the Newtonian center moving along a 
geodesic on S'^ and (the restricted two-body problem). When the curvature 
is small, the pericenter shift is computed using the perturbation theory. We 
also present the results of the numerical analysis based on the analogy with 
the motion of rigid body. 
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1 Problem statement 

We start with the equations of motion for a particle of unit mass on a three- 
dimensional sphere or in a Lobachevsky space (pseudosphere). 

The sphere (pseudosphere L^) can be parameterized using the Cartesian (re- 
dundant) coordinates of the four-dimensional Euclidean space (the Minkovsky 
space M"^) with the constraint 

m = lig.uq'q" ± R') = q) ±R') = o (i) 

where g = diag(l, 1, 1, 1) {g = diag(— 1, 1, 1, 1)) is the corresponding metrics. Here- 
inafter, an upper sign in "±" is used for the sphere and a lower sign is used for the 
pseudosphere. The metrics in M"^ (M"^) generates a metrics in the sphere (the 
Lobachevsky metrics in the pseudosphere L^). 

In terms of the Cartesian coordinates, the Lagrangian for the particle's motion 
in the field of the potential U{q) is 

C = lg,.rr-U{q), 

with the constraint ^ Using the Hamiltonian formalism for systems with constraints 
(Arnold et al. 1993), we get 

^ = UiP^P)-T^)+Uiq). (2) 



2V (q^q) 

rr., , , . dH OH 

Then the canonical equations of motion are q = , p = — . 

op oq 

Two-body problem on (L^)- Consider the two-body problem on the 
curved spaces 5*^ (L^), where bodies are assumed to be point masses. Let these 



masses move in the field of some potential U{qi,q2) (?i,?2 are the coordinates of 
bodies on {L^))- In this particular case the potential energy U depends on the 
distance between two points (this distance is measured along a geodesic). In our 
case, there is a center of mass frame of reference such that the two-body problem can 
be reduced to the problem of the particle's motion in the field of a fixed attracting 
center (i.e. to the Kepler problem in the case of the Newtonian interaction). The 
analogue of the Kepler problem is superintegrable on (L^) (see Kozlov 1994, 
Borisov et al. 1999 and Killing 1885). The generalization of all the Kepler laws 
to spaces of constant curvature is given by Kozlov (1994). But in curved spaces 
the center-of-mass frame of reference does not exist and therefore if the interaction 
between two bodies is Newtonian-like, the two-body problem is not integrable on 
and L3. 

In terms of the Cartesian canonical variables Qa and Pa, the Hamiltonian of the 
two-body system is (the index a denotes the number of the mass m^) 

^=2^ M + 

^ 1 te.P.)few.>-te,<;.P ^^( (3) 

2m2 W2,g2) 

Invariant manifolds. The three dimensional two-body problem is rather com- 
plicated. Therefore by analogy with the planar case, we will examine in detail the 
motion on the invariant submanifolds of the system . The behavior of the system 
on invariant submanifolds allows us to make conclusions about some properties of 
the system (nonintegr ability, stochasticity) in the whole phase space. Nevertheless 



the three-dimensional problem has not been investigated yet. 

The invariant manifolds of the n-body problem in are planes and, similarly, if a 
space is curved, they are spheres S"^ (pseudospheres L^). There is a three-parameter 
family of such manifolds at any point of (L^) (see Borisov et al. 1999). 

Restricted two-body problem. In the Euclidean space a passage to the 
limit is possible in the two-body problem as the mass of attracting center goes 
to infinity while the interaction energy remains finite. There is an inertial frame of 
reference with origin at the "heavier" particle, therefore if the potential is Newtonian, 
the restricted problem is the Kepler problem. 

Consider a similar passage to the limit on (-^^)- In this case the attracting 
center moves along a large circle of the sphere (along a geodesic). The second particle 
(point mass) moves in the field of attracting center and does not affect the motion 
of the attracting center. 

If the origin is at the first particle, the second particle moves in the field of fixed 
center and gyroscopic forces. The Lagrangian of the system is 

^^\^^ + \Y1 ^/'^^m?"^ ^ ^ B^^q^B^pqp - V{q), (4) 

where B = is the angular velocity matrix of the frame of reference. Here, 

B e so(4) (i. e. it is a skew-symmetric matrix) in the case of and B e so(3, 1) in 
the case of L^. 
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2 Bifurcation analysis of the Kepler problem in 
curved spaces 

Let a particle move in the field of Newtonian-like potential on a sphere (pseu- 
dosphere L^). Using spherical coordinates go = RcosO, qi = i?sin6'cos(y9, q2 = 
i? sin ^ sin (y9 cos and q2 = -Rsin^siny^sin-?/', we can write the Hamiltonian as 

Separating variables in (0) gives 

a-ih = Ptb = const, = + . !' = const, (6) 

^ ^ sm (p 



where is the projection of the three-dimensional angular momentum vector M = 
mq X q (here q = (gi, q2, qs)) onto the axis qi, is the squared momentum M"^, E 
is the energy constant. It is easy to see that the vector M is an integral of motion. 
(In the case of the Lobachevsky plane, all the trigonometric functions of 9 should 
be replaced with hyperbolic functions.) 

Let us see how the domain of possible motions on (L^) (hereinafter DPM) 
depends on the energy constant E and the moment constant a^. 

We put r = i?tan6' (r = i?tan6') in (0). If E and are fixed then the DPM 
are defined as follows 



Thus to construct a bifurcation diagram we should consider the quadratic equation 



hr'^ + -fr ^ = 0, (9) 



where h = E ± — The bifurcation set (i.e. the locus of (E,a^) at which the 
2mR^ 

domain of possible motion changes topologically) consists of the curves (see Fig. 1) 



If both roots ri and r2 of are complex (domain I in Fig. 1) the motion is 
impossible. If both roots are real and positive (domain II), the possible values of r 
are given by ri ^ r ^ r2. This implies that a particle moves in the ring 9i ^ 6 ^ O2, 

71 

with < ^^1, < — for 5*^. If the lower root (ri) is negative (domain III), then 
r2 ^ r for the real motion on the Lobachevsky plane and r ^ ri for the motion on 
a sphere (since r is negative if 7r/2 ^ ^ tt). It means that on a body moves 
exterior to the circle 6 ^ 62 and if the space is S"^, a particle moves in the ring 
^1 ^ ^ ^ O2, where < ^1 < 7r/2, 7r/2 < ^2 < vr. 

Note that motions on a sphere are bounded because of compactness S^. Orbiting 
time is always finite. Note also that the "curved" Kepler problems are trajectory 
isomorphic to their plane analogues, as was shown by Serre (see Appell 1891). 

3 Angle-action variables and analogue of Delau- 
nay variables 

Define the action variables in terms of spherical variables 



where the integral is taken over the whole cycle of the period of motion. 

Since = const we have from (fTUj) = = a^. The kinetic energy in terms of 
spherical coordinates on 5*^ is T = -{ped + Pip'P + Piijfp), and on the invariant sphere 

where the particle moves, we have T = -(pe^ + (y<fii^)', here u is a. true anomaly 
(i.e. usual polar angle). Equating these two expressions we get p^pdip = a^du — I^dip- 
The coordinates u and ip change by 27r per one revolution of the orbit. Therefore, 
after integrating we have 

= a^- (11) 

To compute the third integral of (fTn|) . put r = _Rtan6' (r = _Rtan6') and use the 
equation for the orbit r(z/) (see Kozlov 1994; Killing 1885). 



P 



1 + e cos u 



(12) 



where p 
imphes 



where ri 



is a parameter of the orbit, e = W 1 H 5-/^ is the eccentricity. This 



m7 



m7 



' r(l ± r7i?2) 



p p 
r2 



1 + e' 1-e 
We get after integration 

/ 



Ig = V -2mh 



ri 



U2(v/(ri + i?2)(r? + i?2) + i?2 



rir2 



(13) 



for S^, and 



\/ —2mh 



^/iR + n)iR 



r2 



^{R-n){R 



r2 



for L^. 

7 ^ 

Since ri + r2 = — — , rir2 = we get with the exphcit expression of 

h 2mh 

the Hamiltonian 



Similar to the Euchdean space M^, the Hamiltonian depends only on the sum Ig + 

dH 

lio + Iih, i.e. the frequencies Ui = — — , i = 9, ip, tp, corresponding to the variables /e, 

dli 

I^, I^, coincide. This is the case of the complete degeneracy, because all the three- 



dimensional Liouville -Arnold tori foliate into one-dimensional tori i. e. circles. Note 
that unlike the Hamiltonian in the space M.^ (see Markeev 1990), expression 1141 has 
additional terms, which are proportional to — . 

Define new variables L, G, H, I, g, h (analogues of the Poincare variables) 

L = le + + G = I^p + I-ip, H = 
I = (jJe, g = uj^ - ug, h = - lo^. 
In terms of these variables, the Hamiltonian is 



(15) 



m7^ I? , , 



With [TBI and IT^ we have 



"^7 n 



-Eh + ^^2/72 ± l/i?2 ■ 

Equation El implies all the Delaunay variables except / are integrals of motion. 
The angle / is an analogue of the mean anomaly C, and changes uniformly with the 
time I = C = —{t — t). Here r is the time, when the particle passes the pericentre. 



T is the period of orbit revolution, which depends only on the energy constant E 
(see KiUing 1885; Kozlov 1994): 



7 y £;772±i/i?2 

In terms of the angular length of the orbit's major axis a, the energy constant E 



7 7 



L a) 



it!tana V it!tan( 
The Delaunay variables can be expressed in terms of orbit parameters like in the 

planar case as it shown by Markeev (1990) and Demin et al. (1999). Choose the 

angular constants so that if we make gnomonic projection g, h are the images of the 

pericentre parameter and the longitude of the ascending node. Denote them by uj 

and Q,. Let i be the analogue of orbit inclination. This value is equal to the angle 

between the axis q\ and the vector M. 

Express the variables L, G, Z, h in terms of the elements of the orbit p, e, 

T, UJ, VL: 

H = ^/mjp cost, h = Q. 



In the case of the Lobachevsky space L — ^ 777.7 tan (^2) ' 

4 Perihelion shift 

The observation of Mercury's perihelion shift is one of the experiments that proves 
the general relativity theory (GRT) (see Eddington 1963). This shift arises as a 
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result of curving of a space near a gravitating body. Let us prove that in Newtonian 
mechanics, a Keplerian orbit also precesses in a curved space. Although the laws 
of precession in these theories are different. We will take the restricted two-body 
problem as a model problem. This problem is not integrable but if the velocity 
of the heavier particle is low we can analyze the problem using the perturbation 
theory. Here we do not mean to give a new physical justification of the perihe- 
hon shift, already given in GRT and accepted as classical. We just point out that 
some phenomena of the practical Celestial Mechanics admit another interpretations 
(together with the planet nonsphericity, atmosphere refraction and so on). The 
addition of curvature to the classical Newtonian mechanics is an example of such 
interpretations. 

Consider the restricted two-body problem on S"^ {L^)- As usual, the sphere 
(pscudosphere) is assumed to be embedded in IR^ (M^): {q = (x, y, z)\{q, q) = 
-\- y'^ ± — ±E?}. Let an attracting center move along the geodesic on xz plane, 
and we choose the (moving) frame of reference that the attracting center is at the 
north pole of the sphere (pseudosphere) 63 = (0, 0, 1). The Lagrangian of point 
mass (particle of unit mass) is 




here B is the angular velocity matrix of the frame of reference. 











w 



B = 
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Let us assume that the typical size of the domain of motion of the point mass 
is small in comparison with the radius of curvature R. Then we can analyze the 
problem using perturbation theory. Suppose also that the angular velocity of the 
attracting center's motion is small in comparison with the rotation frequency of 
the point mass moving along the corresponding Keplerian orbit. Take the length 
r = i?tan6' (r = i?tan6') and azimuth angle (p as the coordinates on a sphere 
(pseudosphere) and transform UHl as 

if r^0^ \ J w r'^r w"^ r'^ ^ 
£ = - \ o 1 TT H h 2— Tj — coso9=F o- sin oj. (19) 



l^^ipJ i^iP 

V 

Here ^ ~ ^ linear velocity of the motion of the noninertial frame of 

reference. If i? — oo, the problem is reduced to the planar Kepler problem. 

The terms in^J which are linear with respect to the velocity, have the order — 

R^ 

and can not be omitted. To study the evolution of the orbit's shape in the unper- 
turbed Kepler problem we express the equations ^1 in terms of p, u, e, (f. Here, e is 
the eccentricity, u is the longitude of the orbit's pericenter, (p is the azimuth angle, 
p is the orbit's parameter, associated with the energy E of the unperturbed Kepler 
problem by following 

i^ = -^±^. (20) 

The new variables are expressed in terms of coordinates and velocities (hereinafter 7 = 
1) 



p r . esm{(p — u!) r (p 



-r 



1 + e cos(v9 — cu) ' 1 _|_ '"^ \/P ' 11'^^ 

R^ 



Vp- (21) 



Hereinafter, for the sake of simplicity, we don't substitute the expression for r in 
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terms of p, e, u, (p. The Poisson brackets for p, e, uj, (p are 



{e, p} 
and the Hamiltonian is 



Aw pv'^ 

{P, e} = — 2- sin ipsm{(p-uj); 

Aw pr'^ 

{p, cu} = -2 + — 5- sin ip cos(v9 - cj) ; 

K r 

1 - + ^ < 

'^02 Aw pr 

= + ^7 ^ ^' (22) 

{p, ^} = -2v^; 

2 cos(v3 — u) + e + e cos^^ip — u) 



Vp 

sin((/) — uj){2 + ecos{(p — u)) 



-1 9 9 9*2 

1 — e p w r sm (p 



n = -- — -±^±— (23) 

2j9 2i?2 2 r2 ^ ^ 

Expressions 1221 and 1221 imply that when R 00, the variables p, e, uo are slow, p is 

fast. To define the secular change of the orbit's parameters, when R^ r,we neglect 

the terms with order higher than — - and average the equations of motion over the 

R^ 

period of unperturbed motion. Averaging over the period is equal to averaging 

1 f'^'^ 

over ip with a weight function — / f{p)p{p)dp. Here the weight function p is 

2vr Jq ^ 

defined by the derivative p from as p = —. 
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So, we have the system: 

V le-p^l'^ ( 5 e^sincucosct; 

. , P^/^ / , 5 e^sinu;cosa;\ 

f p^/^ /I — 2e^ / 5 

^ = ^^ (i-e^)5/2 ^ + ^ VP ^2 - - cos^ ^ 

The equations 1^ has the integral 

l-e2 



(25) 



2p 

This integral implies that there is no secular change of the energy of the unperturbed 
system EH to this approximation (Laplace's theorem). 

The phase portrait of|^on the surface of the integral EH depends on the param- 

V 

eter h = —=. Figure 2 shows the projection of the trajectories onto {u, e) plane 
V C 



for different values of b. The parameter b describes the ratio of the velocity of the 
attracting center to the characteristic velocity of the particle along the Keplerian 
orbit. The equation implies that the curvature sign determines the direction of 
the motion along the trajectory but not the shape of the trajectory. The value of 
curvature defines the velocity of motion along the trajectory. 

It is clear from the figures that the velocity of perihelion shift depends not only 
on the eccentricity of the orbit but also on the orientation of the orbit with respect 
to the direction of motion of the attracting center. 

If b is small, there exist two stable periodic orbits with the non-zero eccentricity. 
The main axis of the orbits is perpendicular to the direction of the attracting cen- 
ter's motion. At pericenter, the direction of the point mass motion along one of the 
orbit coincides with the direction of the attracting center's motion (^'^ ~ ^) • When 
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the particle moves along another orbit ( ~ ~2~ 



) 



at the pericenter, its direction 



is opposite to the direction of the attracting center's motion. When b increases the 
orbit with ~ ~^ becomes unstable, and if b is sufficiently large, the stable orbit 
with uj = — disappears. 



The projections of the trajectories of the non-averaged system onto {uj, e) plane 
are also shown in the figures. Here, we can see small oscillations (for the variables uj, 
e) near the trajectories of the averaged system (see Fig. 3). 

Remind that standard explanation of the perihelion shift is based on the Schwarzschield 
solution (Eddington 1963) and implies that the shift velocity does not depend on 
the orientation of the orbit u. This is not the case in our problem statement. There 
always exist orbits such that only their eccentricities change their values but their 
pericenters have almost no shift. Moreover, there are fixed points of the system 
()24|) . corresponding to the periodic orbits which do not change their form and ori- 
entation. Note that the Schwarzschield- like metrics can be constructed if boundary 
conditions (at infinity) correspond to the space of constant curvature as it shown by 
Chernikov (1992). And also the restricted two-body problem can be generalized for 
such metrics. It is clear that the velocity of the perihelion shift depends on both u 
and e. 

5 Isomorphism with the spherical top dynamics 

Consider the restricted two-body problem on the sphere S'^ in the general case 
without assumption that the curvature is small. We define the new variables M, 7 
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using the map T*R^ — > e(3), by 

7 = |, M = 7 X p. (26) 

Here, the canonical Poisson brackets {g,, pj} = 6ij are transformed to the Lie - 
Poisson brackets corresponding to e(3) algebra. The equations of motion are 

where the Hamiltonian is 

n=lM' + iM,w) + Uij), f/(7) = -_^L=. (28) 
^ V 7i + 72 

These equations (see Borisov et al. 2001) have two integrals of motion: the area 
integral {M, 7) = C and the geometric integral (7, 7) = 1. In our case C = 0. 
Note that this system describes the motion of spherical top in the potential U{'j) 
and in the field of gyroscopic forces. 

For w = this system is (super)integrable (the Kepler problem on S*^), and in 
this case simple geometric interpretation of the motion exists: the variable M3 = 
const and the projection of the trajectory onto the plane (Mi, M2) is a circle shifted 
from the origin (incidentally, Hamilton noticed a similar thing, in the two-dimensional 
Kepler problem). Indeed in the consequence of {M, 7) = we have M = 7x7 and 

using the equation (fT^ we obtain 

p cos u p sin u 

7i = 72 



a/p^ + (1 + e cos z/)2 a/p^ + (1 + e cos 



1 + e cos z/ p^ez> sin i/ 

73 = =, — — ; rTTTTTTTT = M-i = COnst 

^ v/p2 + (l + ecosz/)2' (p2 + (l + ecosz/)2)3/2 
where i/ is a longitude on the sphere 7^ = 1. Eliminating z> from the last equation 

we can get 

Ml = -p"^M3(e + cos z/), M2 = -p~^M3 sin u 
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For w ^ 0, according to the Liouville- Arnold theorem, an additional integral 
must exist the system to be completely integrable. We will show soon that in the 
general case, the additional integral does not exist (see also Borisov et al. 1999, 
ChernoiVan et al. 1999). 

We construct the Poincare map to study the problem numerically. To construct it 
we use the analogy with the motion of rigid body and choose the Andoyer canonical 
variables (L, G, I, g) according to Borisov et al. (2001) as 

L = Mq, / = arctan — - 
^' VMi 



G = (M, M), g = arccos 



-73 



^1-Mi/{M,M)^ 
Then the Hamiltonian as a function of these coordinates is 



n = -G2 + m;VG2-L2cos/ 



\/G'^ — L? cos g 



2 \/ G"^ sin^ g + L"^ cos^ g 

The equations of motion are canonical: 



j dn . an . an ■ an 

Let us fix the energy level Ti = E and define the Poincare section by the rela- 

TT L 

tion q = —. On this two-dimensional surface we choose the variables — and I as the 
^2 G 

coordinates of the Poincare map (similarly to Borisov et al. 2001). The domain of 

L 



G 



^ 1 and the flow defines 



definition of the variables is compact: / mod 27r, 
corresponding Poincare map. 

By direct substitution into 1221 it is easily proved that 

L(-L, G, g) = -L{L, I, G, g), i{-L, -I, G, g) = -i{L, I, G, g) 

Gi-L, -I, G, g) = G{L, I, G, g), g{-L, G, g) = g{L, I, G, g). 
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So, each trajectory Ci with the initial conditions (Lq, /q; Gq, go) corresponds to a sim- 
ilar trajectory C2 with initial conditions {—Lq, —Iq, Go, go), and each point (L, I, G, g) 
of Ci corresponds to a point {—L, —I, G, g) of C2. This means the Poincare map 
(for the chosen section g = — is central symmetric. The phase portraits for the 
different values of the energy E and parameter w are shown in Fig. 4. 

It is easy to see in the figures that the stochastic layer increases as the energy E. 
This proves that the two-body problem in general case is not integrable. The fixed 
points in Fig. 4 correspond to periodic trajectories of the particle, which play an 
impotent role in the qualitative analysis of the system. 

After the publication of the book by Borisov et al. (1999) and paper by Chernoivan 
et al. (1999), at our suggestion, S. L. Ziglin could prove that the additional mero- 
morphic integral does not exist for the potentials that are the analogous to the 
Newtonian and Hooke interaction for any value of the parameters (see Ziglin 2001 
and Ziglin 2003). 

6 Hill domains and relative equilibrium 

With the equation of motion |2Zl the integral of energy of the restricted problem EHl 
can be written as 

-q^ + U^:{6, ip) = S = const. 



2i?2 

= - sin^ 6^ sin^ V9 - /i cot ^, = ^ > 0, ^ = ^ + 7:, 
2 2 

where 9, ip are the spherical coordinates on 5*^ (it means that 6' is a latitude and if 
is a longitude). 

If the energy E is fixed (therefore, £ is also constant) the domain of motion on 
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S"^ is defined by 

£ -U,{e,ip)^0. (30) 

By analogy with the classical restricted three-body problem (Arnold et al. 1993), we 
will call this domain the Hill domain of the restricted two-body problem on a sphere. 

When the parameters of the system are fixed, the shape of the Hill domains is 
defined by the singularities and critical points of the effective potential U^{9, ip). The 
singularities of U^, are at the poles of a sphere. And, since > and — /i cot 9 > 

— ^0 

oo. Hill domain is always not empty, because near = the inequality (jHUj) holds. 
Each critical point of the function U^{6,ip) corresponds to the equilibrium position 
of the particle (here, the frame of reference rotates with the attracting center). This 
equilibrium position is usually called a relative equilibrium. Note that (we will prove 
this below) in the fixed frame of reference the attracting center moves along the large 
circle and the particle moves along the another circle parallel to the large circle, with 
it being at the same meridian with the attracting center. 
Find the location of critical points by solving the system 



—— = sm cos sm (p -{ 5— = 0, — — = sm t^smc^cosv?. (31) 

09 sin 9 o^p 



We obtain the following results: 

3 TT 

1° if < fi < fi^: = , there are four critical points on the meridians ^ ~ 1^ ^'^'^ 

3 TT 

if = -vr (two points on each meridian) . Their latitudes — < 9i < 92 < n are defined 



by the equation 

lsin2^ + -^ = 0; 
2 sm 9 



2° if fJ^ > fi*, the function Utf{9., ip) has no critical points at all. 



It is easy to see, that in the case 1° the critical points ^— ,^^2^; ^~2~'^^) 

saddle points of the function f/*, and the points (^"^'^i)' ^~2~'^^) strict 

maxima. 

Hill domains for both cases (1° and 2°) are shown in Fig. El It is clear from 
Fig. El that fixed points are in the semisphere opposite to the attracting center. We 
will use linear approximation to investigate the stability of the obtained fixed points 



(relative equilibria). Let the point {(pi,6j), i,j = 1,2, be the corresponding fixed 

vr 3 
2' 2 



point, where (ft G { — ,-7r}, 9j G {6*1,^^2}- According to EH it is very convenient to 



parameterize fi by the latitude of the fixed point: 

/X = — cos 9j sin^ 9j. (32) 

rvr \ TT 
Since /i > for attracting center, 6'j G ^5 Oj) moreover for j = 1 — < 6*1 < 6** and 

these points correspond to the maximum of the effective potential U^, and 6^, < 62 < ir 

to the saddle point. Here, 9^ denotes the value of 9, for which /i EH reaches the 
3^3 

maximum = /i* = . 

16 

Let us introduce canonical impulses pg, corresponding to the spherical angles. 
In the fixed points their values are pg = 0, = ±t(;cos6'j sin6'j. Expand the 
Hamiltonian 1281 in the vicinity of the fixed point up to the second power, using the 
following canonical variables: 

Pe = X, Pip = ±w cos 9j sin 9j + Y, ip = ipi + y, 9 = 9j + x. 
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We obtain 

2 \ sm OjJ \ sm Oj 

+ ]-w^ cos^ Oj ( .^2n + 2/^1 + • • • > ^0 = const, j = 1, 2. 
2 \sm^dj J 

The eigenvalues of the corresponding hnearized system are 



(33) 



. , , 1 — cos 6j — 2 cos^ . , / 1 + cos — 2 cos^ Qj 

1 — COS y 1 + COS 

3 = 1, 2. 

The study of the radical expressions in 1221 gives us the following results: for 
6 = 01 and 6 = 62, Xs^^ are always real and Ai^2 are real for 9 = 61 and purely 
imaginary for 6 = 62- 

This means that 

relative equilibria in the restricted two-body problem on a sphere are always un- 
stable. 

Note that, according to the theorem of central manifold, the existence of two 
purely imaginary eigenvalues for the points (^^5^2^ and ^2^'^^^ results in the 
existence of an unstable (hyperbolic) periodic solution near these points. Fig. U\ 
shows these solutions for different values of energy. 
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Figure 2: Phase portraits of the averaged system 
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Figure 3: Phase portrait of the non-averaged system 
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Figure 4: The Poincare map for the various energy values and for w 



26 



n V 2n 




n V 2n 



Q £ = -0.5 



£ = 



Q £ = 0.525 




-0.2 



^ = 0.5 



^0 ^ = 0.7 



Figure 5: In the case 1° Hill domains (gray shade) for /i = 0.25 < and w = 2.0 
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Figure 6: In the case 2° Hill domains (gray shade) for fi = 0.6 > fi^: and w = 2.0 
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Figure 7: Periodic solutions that appear from the saddle point (7r/2, $2) (a part of 
the sphere is shown in the spherical coordinates). The boundary of the Hill domain 
(see Fig. El for £^ = 0.5) is shown for the critical value of the energy E = when 
other parameters w = 2, 7 = 1. The values of the energy corresponding to shown 
orbits are E=0.0009, 0.01, 0.027, 0.355, 1.0, 1.024 
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